^^Al production from magnetically induced 
extramixing in AGB Stars 
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Abstract 

We discuss nucleosynthesis results obtained following the recent suggestion that 
extramixing phenomena in red giants might be driven by magnetic buoyancy. We 
explore for this model the production of the short-lived radioactive isotope ^^Al 
and of stable light nuclei, considering both the case of the general buoyancy of flux 
tubes and that of the intermittent release of magnetized unstable structures. We 
show that abundant ^^Al can be produced, up to, and above, the highest levels 
measured in presolar grains. This level would be also sufficient to explain the early 
solar system ^^Al as coming from a nearby AGB star of low mass. The case of fast- 
moving instabilities is the most efficient, reaching almost the same effectiveness as 
hot bottom burning (HBB). 
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1 Introduction 



Previous studies have shown that AGB stars are a site where several short- 
lived r adioactivities identified in early solar syste m (ESS) materials are pro- 



duced (IWasserburg et al.l . ll995l : lBusso et al.l . l2003l ). In particular, the possibil- 



ity that one such star polluted the early solar nebula, intro ducing the record 
of nuclei like ^^Al, "^^Ca, ^°Fe, ^°^Pd, ^^^Cs, was explored by IWasserburg et al 



( ll994j ). They found that the polluting star must also experience extra- mixing 
phenor nena, carrying to t he en velope partially H-processed materials rich in 
^^Al (IWasserburg et al.l . 120061 ). These extended mixing episodes were also 
shown to accou nt for the record of C and O isotopes in presolar grains that 
also show ^^Al ( iNoUett et al.l . |2003| ). In the mentioned models circulation was 
treated parametrically, and the physical mechanism was left unspecified. 
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Among the pro cesses invoked to explain extended mixing, rotation played a 



dominant role ( Denissenkov et al. . 20001 ). as it can induce shear instabilities 



and meridian circulation. However, attempts at considering rotational mixing 
i n complete stellar mo dels failed to reproduce the observed abundance changes 
( iPalacios et al.l . l2006l ). This fact promoted a reappraisal of other forms of mat- 
ter transport, like the so-called "tachohaline" mixing in duced by the decrease 
in the molecular weight cre ated by ^He combustion (lEggleton et al 



2008 



Charbonnel and Zahnl . 120071 ). The efficiency of such a mechanism, very im- 



portant in other contexts, seems however to b e in question in the env ironment 



we refer to, i.e. in AGB stars of about 2 Mq (jCantiello et al.l . 120071 ) 



Guided by the existing mode ls of raagneti c buoyancv (jParkerl. Il974 ) and o f 



other magn e tic iri stabilities (iSpruitl . 120011 . 119991 : lEggenberger et al.l . 120051 ) 



Busso et al. ( 20071 ) [hereafter BWNC] recently proposed that matter trans- 
port in red giants be maintained by a magnetic dynamo, creating a system of 
buoyant toroidal flux tubes. In this paper we shall explore in a preliminary 
way some of the nucleosynthesis implications of such an idea. In particular, 
we analyze the synthesis of ^^Al in AGB stars, implementing in a simplified 
way a mass circulation based on BWNC's formalism. 



2 The Model and the Numerical technique 



From the original poloidal field of a rotating star, a toro idal field of s imi 



lar strength can be generated by the dynamo mechanism (jParkerl . Il974j ). In 
physical conditions occurring above the H-burning shell of an evolved star sev- 
eral perturbations can affect a. cylind rical or toroidal flux tube, and they can 



evolve into instabilities (jSpruitl . ll999l ). These are opposed by the magnetic fleld 
aligned with the cylinder, which acts as a stabilizing tension. The complex dy- 
namical behavior of such a tube includes expansion (due to buoyancy), unsta- 
ble oscillations (due to Alfven waves, leading to the development of r2-shaped 



loops ) and torsion due to the Coriolis force (ISpruit and van Ballegooijen 



19821 ). The absence of X-ray emission from surface magnetic activity in AGB 
stars tells us that these magnetized structures, if they do exist, must be dis- 
rupted in the convective layers, thus depositing there their chemical peculiar- 
ities. BWNC showed that efficient buoyancy of magnetized regions requires 
that they achieve very high values of the magnetic fields, up to a few MG. 

It is difficult for a star to maintain very strong fields buried in its interior: if 
the average field is strong, than it is ejected in a short interval of time and a 
whole toroidal fiux tube can be assumed to fioat to the envelope (see Figure 
[H left panel). In this case, maintaining the dynamo for subsequent episodes of 
magnetic activity depends on the possibility that energy is resupplied from the 
convective envelope to the interior and the differential rotation is re-established 
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Fig. 1. Left. A buoyant flux tube, laying in the equatorial plane, with Alfven wave 
oscillations. If the field is strong enough (MG) it floats to the envelope, depositing 
there its composition. Right. Detachment of a magnetized structure due to the 
tilting and reconnection of an unstable 0-shaped loop. The released bubbles move 
rapidly to the envelope, carrying material synthesized close to the H-shell 



( iNordhaus and Balckmaru . l2008l ). If this occurs, than flux tubes emerge, break 
into the envelope and force local matter to be pushed down in a downflow. 
This probably occurs over a wide area, and reaches down to regions where 
p-captures are active. The upward velocity Vu of the buoyancy can be taken 
from BWNC Assuming local mass conservation in each shell crossed by the 
upward and downward motions, the downward velocity would be such that 
Vdfd = Vufu, where fd and /„ are the fractional areas occupied by the upflow 
and by the downflow, respectively. The dimensions of a flux tube as discussed 
by BWNC require that fd/fu be very large (>200) even in the case where the 
activity phenomena related to toroidal fields are constrained within a limited 
range in latitude (in the Sun, this is normally ± 30-50°, corresponding to 
fd = 0.5 — 0.75). In such conditions, H burning has time to occur only in the 
downward movement, due to the lower velocity. When a new magnetic flux 
tube forms, everything is repeated: the net effect is a sort of matter circulation. 
We shall refer to this scheme as Case 1. As discussed in BWNC, the chemical 
changes to be explained require that, in the process, the envelope is mixed 
with a roughly equal mass of partially H-burned material. 



Alternatively, magnetic instabilities might be released intermittently. This cor- 
responds, for example, to the scheme of Figure [1], right panel, in which a kink 
of an Alfven oscillation evolves into an f2-shaped instability of high internal 
fields. It can then be tilted by the Coriolis force, and detached from the rest 
of the flux tube as an independent bubble. We shall refer to this scheme as 
Case 2. We shall present here preliminary results for both cases. 
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3 Effects of magnetic extra-mixing on ^^Al and other isotopes 



We use the model for a star of 2Mq and solar metallicity, previously calcu- 
lated with the FRANEC code using a limited reaction network fo r energy 



production. A mass loss rate as described by the parametrization by iReimers 



(I1975I ). with T] = 0.5 was used. On this basis we recomputed the full nucle- 
osynthesis in the radiative layers above the H shell and followed the mixing of 
processed materials int o the convective env elope. For Case 1, this was done in 



a way similar to that of lNoUett et al.l (120031 ) . but adopting the upward velocity 



given in BWNC and taking = 0.75. For Case 2, we assumed instead that 
every 50 years the phenomenon of Figure [1] (right panel) occurs, carrying to 
the surface H-processed material. The period of 50 years was chosen because 
it is comparable to the magnetic cycles in active stars and to the duration of 



magnetic bursts in evolved stars found by iNordhaud (120071 ). 



In Case 1 we assume a maximum temperat ure for the circulatio n of LogTp 



LogTn — 0.1. This condition was taken from lNoUett et al.l (120031 ) . who showed 
that for higher temperatures excessive energy production would occur. In the 
second case, due to the long duration of the burning at Tp in between two 
episodes of bubble release, we had to reduce Tp not to disturb the stellar 
structure {ALogTp = 0.15). In several runs we explored the range of mass 
circulation fr o m M =10~'' to M=10~^ MQ/yr, following the indications by 
NoUett et all (l2003h that this is the interval in which the isotopic shifts ob- 



served in presolar grains of AGB origin can be explained. We notice that, 
although the value of M is not self-consistently obtained, but rather adopted 
from the quoted work, there are specific relations linking the mass circulation 
to the magnetic field B (at least in Case 1). Indeed, for flux conservation, call- 
ing a the radius of the flux tube, M oc Ba^. On the other hand, the buoyancy 



velocity v is related to the field by the relation v oc B (see iBusso et al.l . l2007l . 
Eq.4), therefore M oc B x a2 . 

Tube dimensions were assumed from the solar case; then increasing B is re- 
quired to obtain higher M values. Similarly, lower mass circulations would 
imply lower fields, as indeed found by BWNC for RGB stars. While small 
vales of the mass circulation are certainly possible, it seems unlikely that a 
much faster transport (say, by a factor 100) can be magnetically supported, 
unless the tube dimensions are very different from the solar case. Indeed, the 
field values obtained by BWNC for M ~ 10~^ M^/yr are already quite high. 
In the scheme of Case 2, the analogy with solar magnetic tubes is less im- 
portant, because the achieved envelope abundances are essentially controlled 
only by the long burning periods (tens of years) at Tp. In this case, magnetic 
buoyancy simply transports in discrete bubbles matter synthesized by the H 
shell, but continuously resupplying the fuel with material descending from 
the envelope (by mass conservation). High B values are required only in the 
buoyant bubbles, not in the buried magnetic structures, and a faster mass 
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transport is a priori possible. 



In general, computing extra-mixing requires describing nuclear physics phe- 
nomena coupled with dynamics. This can be obtained by solving a system of 
differential equation of the type: 



dNi _ dNi dNi dM OR 
"df ~~df^ dM dR 'dt 



where the partial time derivative due to nucleosynthesis is: 
dN- 

-NpNi\i,p + ArpA/:,_iAi_i,p - NiX^ + NAd (2) 



dt 



Here the par ameters A are t he rea ction rates, which we take from the NACRE 
compilation (lAngulo et al.l . Il999l ). In practical cases, the phenomenon might 
become simpler than that. For example, in the fast upward flux the time scale 
for motion is smaller than the one for burning, and no nucleosynthesis occurs. 
On the contrary, in a sufficiently slow descending flux everything depends only 
on the path integral o f reac tion rates, as in previous extramixing calculations 
by IWasserburg et al.l (Il995l ) 



In our work we first recomputed the " static" H burning through the radiative 
region for the whole AGB phase, with an extended network going from H to 
^^S. Then we applied separately our mixing schemes for Cases 1 and 2. 

Figure [2] shows preliminary results thus obtained, in terms of abundances 
of affected isotopes in the envelope, as a function of the time spent on the 
thermally-pulsing AGB stage. The changes in abundances introduced by the 
third dredge-up have for the moment been neglected. 

Preliminary conclusions drawn from Figure [2] can be summarized as follows: 

• Case 1 (Figure [21 left panel) gives results very similar to what was previously 
obtained in parameterized extramixing (or c ool bottom processin g, CBP, 
as this phenomenon is sometimes called) by iNoUett et al.l (120031 ). ^^Al is 
produced at a moderate efficiency, reaching final abundance ratios in the 
envelope in the range ^^Al/^^Al = 10~^ to 10~^ ( depending on the choice o f 
parameters). This corresponds to values found by IWasserburg et al.l (120061 ). 
^^O is virtually untouched (i.e. we have only CN cycling), while heavier 
oxygen isotopes are largely modified, covering the region where presolar 
grain data lay. 

• Case 2 (Figure [21 right panel) shows a somewhat different burning process. 
The abundances of nuclei modified at high T (like ^^O) are more affected, 
and ^^Al/^^Al reaches values close to 0.1. By contrast, fragile nuclei like ^^O 
and ^^N, already burning at low T, are less affected than in Case 1, where 
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Fig. 2. Left: envelope abundance evolution along the TP-AGB phase for nuclei 
affected by CBP, in a star of IMq and Zq^ in Case 1 (circulation). Right: envelope 
abundances in Case 2 (detached bubbles) . Here ti is the the time spent between the 
release of two successive bubbles 

they have more time to be consumed. The high abundance of ^^Al found 
makes this case interesting for the possibility that the encounter of an AGB 
star be at the origin of the early solar system concentration of ^^Al. 
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